A numerical study of squeeze-film damping in MEMS-based structures including rarefaction effects

In a variety of MEMS applications, the thin film of fluid responsible of squeeze-film damping results to be rarefied and, thus, not suitable to be modeled though the classical Navier-Stokes equation. The simplest way to consider fluid rarefaction is the introduction of a slight modification into its ordinary formulation, by substituting the standard fluid viscosity with an effective viscosity term. In the present paper, some squeeze-film damping problems of both parallel and torsion plates at decreasing pressure are studied by numerical solving a full 3D Navier-Stokes equation, where the effective viscosity is computed according to proper expressions already included in the literature. Furthermore, the same expressions for the effective viscosity are implemented within known analytical models, still derived from the Navier-Stokes equation. In all the considered cases, the numerical results are shown to be very promising, providing comparable or even better agreement with the experimental data than the corresponding analytical results, even at low air pressure. Thus, unlike what is usually agreed in the literature, the effective viscosity approach can be efficiently applied at low pressure regimes, especially when this is combined with a finite element analysis (FEA). SOMMARIO. In molte applicazioni MEMS, il sottile strato di fluido responsabile del fenomeno dello squeezefilm damping risulta essere rarefatto e, quindi, non modellabile mediante l’equazione classica di Navier-Stokes. Il modo più semplice per tener conto della rarefattezza del fluido consiste nell’introduzione di una piccola modifica all’interno della sua formulazione ordinaria, ovvero nella sostituzione della viscosità del fluido con una viscosità effettiva. Nel presente lavoro, vengono presi in considerazione alcuni problemi di squeeze-film damping, riguardanti casi in cui le piastre coinvolte sono dotate di moto normale traslatorio e casi in cui esse sono dotate, invece, di moto torsionale. Tali problemi vengono risolti mediante una modellazione numerica 3D dell’equazione di Navier-Stokes, in cui la viscosità effettiva viene calcolata mediante apposite espressioni già note in letterature. Inoltre, queste stesse espressioni sono utilizzate all’interno di modelli analitici già conosciuti, derivati anch’essi dall’equazione di Navier-Stokes. In tutti i casi qui considerati, i risultati numerici sono molto promettenti, in quanto essi sono caratterizzati da uno scostamento dai dati sperimentali uguale o inferiore a quello fornito dai corrispettivi risultati analitici, anche a basse pressioni. Quindi, nonostante quanto sia usualmente asserito in letteratura, il metodo della viscosità effettiva può essere efficientemente adottato anche a bassi regimi di pressione, specie se in combinazione con analisi numeriche agli elementi finiti (FEA).